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Abstract 

We present an analysis of Galactic halo structure and substructure 
traced by ~ 500 RR Lyrae stars in the SDSS stripe 82 region. The main 
result of this study is a 2D map of the Galactic halo that reaches distances 
of 100 kpc and traces previously known and new halo substructures, such 
as the Sagittarius and Pisces tidal streams. We also present strong direct 
evidence, based on both RR Lyrae and main sequence stars, that the halo 
stellar number density profile significantly steepens beyond 30 kpc from 
the Galactic center. Using a novel photometric metallicity method that 
simultaneously combines data for RR Lyrae and main sequence stars, 
we show that the median metallicity of the Sagittarius trailing stream 
in the SDSS stripe 82 region is [Fe/H] = —1.2 ± 0.1 dex. In addition to 
these results, we will present ongoing work on a 3D map of the Galactic 
halo constructed using a sample of ~ 4000 RR Lyrae stars that covers 
8500 deg 2 of northern sky and probes up to 30 kpc from the Sun. 

1. Galactic Halo as the Rosetta Stone for Galaxy Formation 

Studies of the Galactic halo can help constrain the formation history of the Milky 
Way and the galaxy formation process in general. For example, state-of-the-art sim- 
ulations of galaxy formation predict numerous substructures, s uch as tidal tails and 



streams, in halos of Milky Way -size galaxies (e.g., see Fig. 14 in I Johnston et all 12008 



ay-; 

or Fig. 6 in lCooper et al.ll2010l ). The amount, morphology, kinematics, and chemical 



compos ition of these substructures depend on th e accretion history of the simulated 



galaxy (1 Johnston et al.ll2008l ; ICooper et al.ll2010l ). Therefore, if we map the substruc- 



tures in the Galactic halo and compare the resultant maps with simulations, we will 
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be able to constrain the formation history of the Milky Way. 

Maps of the Galactic halo that probe galactocentric distances (R gc) of 



20-30 



kpc have alread y been made using near turn-off main sequence stars (jJuric et a. 



Bell et al.ll2008l ) selected from the Sloan Digital Sky Survey (SDSS; lYork et al. 



2008, 
2000h . 



Unfortunately, the largest discrepancy between simulated halos with different formation 
histories occurs at distances beyond 30 kpc (e.g., compare panels in lCooper et al.ll2010 



Fig. 6). Therefore, to discriminate between different formation histories we need to 
map the Galactic halo beyond 30 kpc and to do that we need to use tracers that are 
brighter than main-sequence turnoff stars. 



2. Mapping the Galactic Halo with SDSS Stripe 82 RR Lyrae 

Stars 

To map the Galactic halo beyond 30 kpc, we used RR Lyrae stars selected from 
the SDSS stripe 82 region (20 h 32 m < a J20 oo.o < 04 h 00 m , -1.26° < 5 J20 oo.o < +1-26°, 
~ 280 deg 2 ). This region has been repeatedly observed during SDSS-I and SDSS-II 
parts of the survey and on average, there are about 30 observations per source. The 
photometric precision of SDSS Stripe 82 dat a is ~ 0.02 mag at the bright end (r < 18) 
and ~ 0.05 mag at r ~ 21 (ISesar et al.ll2007l ). For comparison, the apparent magnitude 
of an RR Lyrae star at 100 kpc is r ~ 21. 

The candidate RR Lyrae stars were pre-selected using SDSS colors and using some 
low-level variability statistics, such as rm s scatter and chi-sq uare per degree of freedom 
(for details on the selection algorithm see Sesar et al. 2010al ). The period of variability 
was found for ea ch candidate using an implementation of the Supersmoother algorithm 
f Reimannl 1994 ). the candidate light curve was folded using that period and a custom 
set of ugriz template light curves was fitted to the period-folded light curve for identi- 
fication. In total, we have identified about 400 type-a& (RRab) and about 100 type-c 
(RRc) RR Lyrae stars in the SDSS stripe 82 region. 



To calculate distances for RRab stars we used the IChaboyerl (119991 ) My — [Fe/H] 
relation 

M v = (0.23 ± QM)[Fe/H\ + (0.93 ± 0.12) (1) 
and assumed [Fe/H] = —1.5 dex (median halo metallicity; llvezic et al.ll2008l ) as the 



metallicity of all RR Lyrae stars in our sample 1 . This relation gives My = 0.6 for the 
absolute magnitude of RRab stars in the Johnson V band. The resulting fractional un- 
certainty in distance due to unknown metallicity, evolutionary effects, and photometry 
is ~ 5%. RRc stars were not used in the rest of this study. 

A Bayesian number density estimator (llvezic et al.ll2005l ) was applied on the spatial 
distribution of RR Lyrae stars and the calculated number densities of RR Lyrae stars 
were compared to values predicted by a smooth, oblate halo model with a density 



1 Spectroscopic metallicities were not available for the majority of stars in our sample. 
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Figure 1. — Halo substructures traced by RR Lyrae stars in SDSS stripe 82. The heliocentric 
distance in kpc is in the radial direction and angle indicates equatorial right ascension. The red 
color indicates substructures, green color shows the smooth halo, and the blue color shows regions 
where there are almost no stars. The Pisces Overdensity /Strea m is located ~ 80 k pc from the Sun at 
R.A. ~ 355 deg, the trailing arm of the Sagittarius tidal stream (jlvezic et al.ll2003l ) is passing through 
the stripe 82 plane at (R.A., helio. dist.)=(~ 30 deg , 25 kpc), a nd the Hercules- Aquila Cloud is 
located at (~ 330 deg, < 25 kpc). Figure adapted from lSesar et all (|2010aD Fig. 11. 



power-law slope of -2.7 (best-fit halo model from lJuric et all 120081 ). Fig. [T] shows 
the comparison of observed and model-predicted number densities on a logarithmic 
scale. The green regions are in agreement with the smooth halo model, red regions are 
overdense by a factor of 10, and blue regions are underdense by a factor of 10 compared 
to the model. 

Within about 30 kpc, the observed halo follows the oblate power-law model. Two 
substructures are clear ly visible: the Sgr trailing arm and the Hercules-Aquila Cloud 
(IBelokurov et all 120071 ). Beyond 30 kpc, the model predicts more stars than what is 
actually observed, indicating that the observed number density profile steepens. This 
steepening is also pr e sent w hen main-sequence stars are used as tracers, as shown in 
Sesar. Juric &: Iveziq (120111 ) Figs. 9 and 10 and in Fig. [2] below. 



The overdensity seen in Fig. [T] at about 80 kpc was first reported in ISesar et al 



f|2007h ("J" clump) and was l ater r enamed the Pisces overdensity by IWatkins et al. 



(120091 ) . iKollmeier et all (120091 ) and ISesar et all (l2010bf ) have observed this substruc- 
ture spectroscopically and have detected two velocity peaks, suggesting that the sub- 
structure may be a tidally di srupted dwarf galax y. A more extended view of this 
substructure was provided by ISharma et al. (bold ) who used M-giants selected from 
2MASS. 
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Figure 2. — Halo substructures traced by near turnoff main-sequence stars in SDSS stripe 82. 
The distance limit in this map is at 40 kpc from the Sun. Note the Sagittarius tidal stream, Hercules- 
Aquila Cloud, and the steepening of the halo density profile beyond ~ 30 kpc. Figure adapted from 
iSesar et al.l (|2010al ) Fig. 24. 



2.1. A New Photometric Metallicity Method 

As shown in Figs. [I] and El we have detected the Sagittarius dSph tidal stream 
(trailing arm) in SDSS stripe 82 as an overdensity of RR Lyrae and main-sequence stars. 
These detections have allowed us to develop and test a new photometric metallicity 
method that can be used when SDSS u band observations are not available. 

As shown in Fig. El for fixed apparent magnitudes, the distance modulus of RR 
Lyrae and main-sequence stars changes with metallicity, but with opposite signs. For 
RR Lyrae stars, the distance modulus decreases with increasing metallicity, as dot- 
dashed lines show. The distance modulus of main-sequence stars, on the other hand, 
increases with increasing metallicity, as dashed lines show. If we assume that a clump of 
main-sequence stars belongs to the same substructure as the clump of RR Lyrae stars, 
we can solve for metallicity and distance modulus. In the case of Sagittarius trailing 
arm this method gives median metallicity of a bout -1.2 dex , and agrees well with 
the median spectroscopic metallicity obtained by lCarlin et al.l ( 120101 ) ([Fe/H] = —1.15 
dex). 

The u—g color of main-sequen ce stars can also be used as an estimator of metallicity, 



as detailed in llvezic et al.l (120081 ) . Using this method and main-sequence stars in SDSS 
stripe 82, we estimated —1.2 dex as the metallicity of the Sagittarius trailing arm. 
However, the Ivezic et al. photometric metallicity method requires u — g observations 
of main-sequence stars. These observations are difficult to obtain and some current 



Mapping the Galactic Halo with RR Lyrae Stars 



17.8 
1 7.7 
17.6 
17.5 
1 7.4 
^ 17.3 
2 17.2 

Q 

1 7.1 
17 
16.9 
16.8 



Sgr tidal stream [Fe/H] and DM in Stripe 82 



1 1 1 1 r 
RR Lyrae 



n 1 1 r 



1 ' 1 1 "--7 
/ 

/ -I 

/ 

/ 

/ 

/ _/ 
/ /- 
/ / 
/ / 
/ / 
/ / 
/ / 
' / 



MS 



' i ^ 

A 



-A- 



/ 



u-g 



j i i <d_ 



J I L 



-2 



-1.5 -1 
[Fe/H] 



-0.5 



Figure 3. — A summary of the constraints on the distance and metallicity of the Sagittarius 
dSph tidal stream (trailing arm). The dot-dashed lines show a constraint obtained from the mode of 
the apparent magnitude distribution of RR Lyrae stars, with ±0.1 mag adopted as the uncertainty. 
Diagonal short-dashed lines are analogous constraints obtained from the median apparent magnitude 
of main-sequence stars with 0.4 < g — i < 0.5. The vertical lon gdashed line s mar k the median 
photometric metallicity for main sequence stars obtained using the llvezic et all ()2008l ) photometric 
metallicity method, with ±0.1 dex adopted as the uncertainty. All three constraints agree if the 
distance modulus is DM = 17.2 mag and [Fe/H] = —1.2 dex. 

and upcoming sky surveys will not have w-band detections at all. The advantage of 
the above method is that it does not need w-band observations, and therefore it can be 
used by surveys such as PTF, PanSTARRS and DES that will not have observations 
in the u band. 
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3. Mapping the Galactic Halo with LINEAR and PTF RR 

Lyrae Stars 

The biggest disadvantage of SDSS stripe 82 data is its small sky coverage (only 
~ 1% of the sky), meaning that some of the results and conclusions presented here 
may not be representative of the entire halo. To address this issue, we have selected 
about 4000 RRab stars from the LINEAR survey (covering 10,000 deg 2 of sky) and 
have started to map the Galactic halo in 3D in order to verify the above results based 
on SDSS stripe 82 data with much smaller sky coverage (~ 300 deg 2 ). 

LINEAR stands for the Lincoln Near-Earth Asteroid Research. This ongoing survey 



(jStokes et al.ll2000f ) has been observing the northern sky since 1998 and is a premier 
source for variability studies. Unfortunately, detections of asteroids were more impor- 
tant to this survey than photometry so the quality of original data is not as good as 



i n SD SS stripe 82. Following a procedure similar to the one described in ISesar et al. 



(120061 ) . we have recalibrated LINEAR photometry using SDSS data and have created 
a database that contains 8.2 million objects with more than 200 observations per ob- 
ject. A sample of about 4000 RRab s t ars has been selected from this dataset using the 



selection algorithm from ISesar et al.l ( ]2010af ). This sample is 90% complete up to 30 
kpc, covers 10,000 deg 2 of sky, and the contamination (fraction of non-RR Lyrae stars) 
is less than 2%. 

In order to detect substructures in the region of halo pro bed by LINEAR RR Lyrae 



stars, we have applied the Enlink group-finding algorithm ( ISharma fc Johnston! 120091 ) 
to our sample of LINEAR RRab stars. The Enlink algorithm has identifed 6 significant 
halo groups based only on positions of RR Lyrae stars, and we show these groups in 
Fig. |5j Even though most of these groups can be associated with g lobular clusters or 



known halo streams (e.g., Virgo Stellar Stream; iDuffau et al.l 120061 ) . the morphology 
of some groups seems to indicate that they are tidal streams, possibly originating from 
some of the globular clusters. 

Even though the LINEAR survey covers a much wider area of the sky than SDSS 
stripe 82, it is quite shallow (r < 18) and RR Lyrae selected from it probe only the 
inner part of the Galactic halo (within 30 kpc from the Sun). What is needed is a 
multi-epoch survey that covers a wide area of the sky and is deep at the same time. 
The survey that satisfies these criteria and is the Palomar Transient Factory. 

PTF is a wide-area, two-band (SDSS-g' and Mould-R filters), deep (R ~ 21, single- 
epoch; R ~ 23, co-added) survey aimed at systematic exploration of the optical tran- 
sient sky. It will provide accurate photometry (with systematic uncertainties < 0.02 
mag) for more than 10,000 deg 2 of sky in two cadences: i) 5DC — two 60s exposures 
separated by one hour and repeated after 5 days, and ii ) DyC — a dynamical cadence 
designed to explore transient phenomena on timescales longer than 1 min and shorter 
than 5 days. These are optimal cadences for the observation of RR Lyrae stars (period 
of pulsation ~ 0.6 days) and will ensure a dense phase coverage of light curves. This 
survey is ongoing and will provide 60 epochs by the end of the survey, enough for an 
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Figure 4. — The median number of observations per object in the recalibrated LINEAR dataset 
as a function of equatorial J2000.0 right ascension and declination coordinates. The values are color- 
coded according to the legend, with values outside the range saturating. The dashed lines show ±10° 
of the ecliptic plane (solid line) and the Galactic plane is shown as a dotted line. The average number 
of observations per object within ±10° of the ecliptic plane is ~ 460, and ~ 200 elsewhere. 



efficient selection of RR Lyrae stars. 



4. Conclusions 



We have presented the most complete sample of RR Lyrae stars identified so far in 
the SDSS stripe 82 data set, consisting of 379 RRab and 104 RRc stars. Our visual 
inspection of single-band and color light curves insures that the sample contamination 
is essentially negligible. Although th e sky area i s relatively small co mpared to other 



recent surveys, such as iKeller et al.l ( 120081 ) and iMiceli et al.l (120081 ). this RR Lyrae 



sample has the largest distance limit to date (~100 kpc). 

The high level of completeness and low contamination of our resulting sample, as 
well as the more precise distance estimates, enabled a more robust study of halo sub- 
structure than was possible in our first study. Our main result remains: the spatial 
distribution of halo RR Lyrae at galactocentric distances 5-100 kpc is highly inho- 
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Figure 5. — A 3D rendering of halo substructures detected using LINEAR RRab stars. The 
points show RRab stars clustered by Enlink into 6 significant groups (there are two red, blue, and 
white clusters). The Sun's position in the Galactic disk is indicated by a yellow sphere. The green 
spheres show positions of s ome globular cluster s and the blue sphere shows the position of the Virgo 
Stellar Stream according to lDuffau et al.l ()2006f ) . The spatial extent of some groups suggests that they 
may be tidal streams. 



mogeneous. A comparison of the observed sp atial distribution of RR Lyrae stars and 
main sequence stars to the lJuric et al. (120081 ) model, which was constrained by main 
sequence stars at distances up to 20 kpc, strongly suggests that the halo stellar number 
density profile steepens beyond ~30 kpc. While various indirect evidence for this be- 
havior, bas ed on kinematics of field stars, globular clusters, and other tracers has been 
published f Carollo et al. 20071 : and references therein), our samples provide a direct 
measurement of the stellar halo spatial profile beyond the galactocentric distance limit 
of ~30 kpc. A si milar steepening of the spatial profile was detected using candidate 



RR Lyrae stars bv Keller et al.l ( 120081 ) . 



We introduced a novel method for estimating metallicity that is based on the ab- 
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solute magnitude vs. metallicity relations for RR Lyrae stars and main sequence stars 
(calibrated using globular clusters). This method does not require the w-band photom- 
etry, and will be useful to estimate metallicity of spatially coherent structures that may 
be discovered by the Dark Energy and Pan-STARRS surveys. While the existing SDSS 
data are too shallow to apply this method to the Pisces stream, we used it to obtain 
a metallicity estimate for the Sgr tidal stream that is consistent with an independent 
estimate based on the photometric w-band method for main sequence stars. Our result, 
[Fe/H] = —1.2, with an uncertainty of ~0.1 dex, strongly rules out the hypothesis that 
the trailing arm of the Sgr dSph tidal st ream has the s a me m et allicity as ha l o field 
stars ([Fe/H] = -1.5), as suggested bv IWatkins et all ( 120091 ). lYannv et all feOQOh 
detected peaks at [Fe/H] = —1.3 and at [Fe/H] = —1.6 using SDSS spectroscopic 
metallicities for blue horizontal-branch (BHB) stars from the trailing arm. However, 
both the Yanny et al. and Watkins et al. results could be affected by systematic er- 
rors in the SDSS metallicity scale for BHB stars. Our results suggest that the SDSS 
metallicity scale for BHB stars could be biased low by about 0.3 dex (relative to the 
SDSS metallicity scale for main sequence stars, and assuming that RR Lyrae and main 
sequence stars from the Sgr tidal stream in stripe 82 area have the same metallicity 
distributions). 

Simulations by Bullock &: Johnston! ( 120051 ) and IJohnston et al.l (120081 ) predict that 
there should be a difference in the chemical composition between stars in the inner 
halo that was built from accretion of massive satellites about 9 Gyr ago, and outer 
halo dominated by s tars coming from dSph sa tellites that were accreted in the last 
5 Gyr (see Fig. 11 in iBullock Sz Johnston l2005h . These accreted dSph satellites were 
presumably more metal-poor than the massive satellites accreted in earlier epochs 



Robertson et al.l 120051). Other recent simulation studies support these 



De Lucia fc Helmil ( 120081 ) find evidence in their simulation 



(see Fig. 3 in 
conclusions. For example, 
for a strong concentration of (relatively) higher metallicity stars at distances close to 
the Galactic center, and the pres ence of (relatively) lo wer metallicity stars at distances 
beyond 20 kpc from the center. IZolotov et al.l ( 120091 ) find from their simulations that 
their inner halos include stars from both in-situ formed stars and accreted populations, 
while their outer halos appear to originate through pure accretion and disruption of 
satellites. 

Simulations indicate that high surface brightness substructures in the halo originate 
from single satellites, typically m assive dSph which tend to be accreted over the last few 
Gyr (IBullock fc Johnston! 120051 ). and these massive ga laxies are expected to be more 



more metal-rich than halo field stars ( iFont et al.l 120081 ) . The results from llvezic et al. 



(120081 ) and the results presented here seem to support this prediction. The inner halo 
has a median metallicity of [Fe/H] = —1.5, while at least two strong overdensities have 
higher metallicities - the Monoceros stream has [Fe/H] = —1.0, and for the trailing 
part of the Sgr tidal stream we find [Fe/H] = —1.2. We emphasize that these three 
measurements are obtained using the same method/calibration and the same data set, 
and thus the measurements of relative differences are expected to be robust. 
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Our result that (inner-) halo stellar number density profile steepens beyond ~30 kpc 
is limited by the relatively small distance limit for main sequence stars (35 kpc), the 
sparseness of the RR Lyrae sample (~500 objects), and the small survey area (~300 
deg 2 ). Ideally, the halo stellar number density profile should be studied using numerous 
main sequence stars detected over a large fraction of sky. To do so to a distance limit 
of 100 kpc, imaging in at least g and r bands (or their equivalent) to a depth several 
magnitudes fainter than the co-added SDSS stripe 82 data is required (r > 25). Pan- 
STARRS, the Dark Energy Survey and LSST are planning to obtain such data over 
large areas of sky. The LSST, with its deep w-band data, will also exten d metallicity 



mapp ing of field main sequence stars over half of the sky in the south; see llvezic et al. 
(120081 ) for details. For substructures to be potentially discovered in the north by Pan- 
STARRS, the method presented here can be used to estimate the metallicity of spatially 
coherent structures even without the w-band data. 
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